Detection of a persistent-current qubit by resonant activation 
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We present the implementation of a new scheme to detect the quantum state of a persistent- 
current qubit. It relies on the dependency of the measuring Superconducting Quantum Interference 
Device (SQUID) plasma frequency on the qubit state, which we detect by resonant activation. With 
a measurement pulse of only 5ns, we observed Rabi oscillations with high visibility (65%). 
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Superconducting circuits containing Josephson junc- 
tions are promising candidates to be the building blocks 
(qubits) for future quantum computers Q . Coherent dy- 
namics has been observed with several different qubit 
designs. Complex sequences of rotations on the Bloch 
sphere could be performed in NMR-like experiments : 
Rabi, Ramsey, spin-echo, composite pulse sequences were 
demonstrated @, 0, 0.HJ. I* 1 these experiments, the de- 
tection of the qubit energy state is done after the coher- 
ent operations have been accomplished by reading out a 
macroscopic detector with an output that is correlated 
to the qubit energy eigenstate. By averaging the sig- 
nal delivered, one infers the qubit eigenstate occupation 
probability. 

A hysteretic DC SQUID has been succesfully em- 
ployed to measure the state of a frux-qubit @ (a simi- 
lar scheme was used in Q in the case of a split Cooper 
pair box). In this experiment, a fast dc current pulse 
(DCP) sent through the SQUID induced switching to 
the finite-voltage state conditional on the qubit state. 
In this Rapid Communication, we propose and demon- 
strate a new read-out method, faster than the previous 
one, which does not require a large bandwidth line to 
inject current through the SQUID and thus will allow a 
stronger filtering of the bias lines. The method is based 
on the dependence of the SQUID plasma frequency on 
the qubit 's state measured by a resonant activation mi- 
crowave pulse (RAP). We compare Rabi oscillations ob- 
tained by the two schemes in exactly the same conditions, 
and show that the RAP method significantly increases 
the contrast. 

The persistent-current qubit consists of a micron-size 
superconducting loop intersected with three Josephson 
junctions |6j. It is threaded by a magnetic flux <f> x gen- 
erated by an external coil. When the total phase across 
the three junctions j q is close to tt (meaning that <$> x is 
close to $o/2), the loop has two low-energy eigenstates 
(ground state and excited state 1) well separated from 
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the higher-energy ones, and can thus be used as a qubit. 
They are linear combinations of two states | f ) and | |) 
which carry an average persistent-current of +I p and — I p 
and are coupled by tunnelling. The effective hamilto- 
nian reads H = —h(e<r z /2 — Aa x /2) in the {| f),| J.)} 
basis, where hA/2 is the tunnelling matrix element be- 
tween the two basis states, and e = 2I p (Q x — <I>o/2)//i. 
This leads to an energy separation between the two states 
Ex-E = hf q = hVA 2 + e 2 0. 

The average current in state i [i = 0, 1) 
Jj = — < dH/d$ x >i can be computed from the previ- 
ous equations. It is shown in figure^ for the parameters 
of our sample. The qubit is inductively coupled to a hys- 
teretic DC-SQUID whose critical current depends on li 
by an amount MliCe, where M is the mutual inductance 
between the qubit and the SQUID, and a = die /d§ x is 
the slope of the SQUID modulation curve. The require- 
ment imposed upon any detection scheme is to detect the 
small (about 2%) variation in the SQUID critical current 
associated with a transition between the qubit states in 
a time shorter than the qubit 's energy relaxation time 
T\. A first method was demonstrated in @- It con- 
sists of applying a short DCP to the SQUID at a value 
lb during a time At , so that the SQUID will switch 
out of its zero- voltage state (0 ; @) with a probability 
Psw{h)- For well-chosen parameters, the detection effi- 
ciency D = max h \P} w - P°J, where P° w (P} w ) is the 
switching probability if the qubit is in the ground (ex- 
cited) state, can approach 1. The switching probability 
then directly measures the qubit 's energy level popula- 
tion. 

The Al/AlO x /Al tunnel junctions as well as the 
SQUID and qubit's aluminium loops were obtained by 
e-beam lithography and shadow-evaporation on an oxi- 
dized silicon substrate. The SQUID was shunted by an 
on-chip aluminium capacitor and connected to the cur- 
rent injection line through a 400 fl on-chip gold resistor 
(see figure^,). Both current and voltage line were fil- 
tered by copper-powder filters The voltage across 
the SQUID was amplified at room-temperature and sam- 
pled by an acquisition card. By repeating the measure- 
ments typically 10000 times at a rate of 25 kHz we derive 
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FIG. 1: (a) AFM picture of the SQUID and qubit loop (the 
scale bar indicates 1 (im). Two layers of Aluminium were 
evaporated under ± 20 ° with an oxidation step in between. 
The Josephson junctions are formed at the overlap areas be- 
tween the two images. The SQUID is shunted by a capacitor 
Csh = 12 pF connected by Aluminium leads of inductance 
L = 170 pH (solid black line). The current is injected through 
a resistor (grey line) of 400 SI. (b) DCP measurement method 
: the microwave pulse induces the designed Bloch sphere ro- 
tations. It is followed by a current pulse of duration 20 ns, 
whose amplitude h is optimized for the best detection effi- 
ciency. A 400 ns lower-current plateau follows the DCP and 
keeps the SQUID in the running-state to facilitate the voltage 
pulse detection, (c) Larmor frequency of the qubit and (in- 
sert) persistent-current versus external flux. The squares and 
(insert) the circles are experimental data. The solid lines are 
numerical adjustments giving the tunnelling matrix element 
A, the persistent-current I p and the mutual inductance M. 

the switching probability. Pulsed microwave irradiation 
was applied to the qubit through a coaxial cable attenu- 
ated by 20 dB at T = 1.5 K. The sample was enclosed in 
a copper box thermally anchored to the mixing chamber 
of our dilution refrigerator (base temperature 25 mK) to 
shield it from RF noise. 

The parameters of our qubit were determined by fit- 
ting spectroscopic measurements with the above formu- 
lae. For A = 5.855 GHz, I p = 272 nA, the agreement is 
excellent (see figure^). We also determined the coupling 
constant between the SQUID and the qubit by fitting 
the qubit "step" appearing in the SQUID's modulation 
curve (see insert of figure HJ;) and found M — 20 pH. 
We first performed Rabi oscillation experiments with the 
DCP detection method (figure^)). We chose a bias point 
& x , tuned the microwave frequency to the qubit reso- 
nance and measured the switching probability as a func- 
tion of the microwave pulse duration r mw . The observed 
oscillatory behavior (figure^) is a proof of the coherent 
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FIG. 2: (a) Rabi oscillations at a Larmor frequency f q — 
7.15 GHz (b) Switching probability as a function of current 
pulse amplitude I without (closed circles, curve Pg W (I)) and 
with (open circles, curve PJL (-/&)) a n pulse applied. The solid 
black line P° h (I b ) is a numerical adjustment to Ps W (Ib) as- 
suming escape in the thermal regime. The dotted line (curve 
Pth(h)) is calculated with the same parameters for a critical 
current WOnA smaller, which would be the case if state 1 was 
occupied with probability unity. The grey solid line is the 
sum0.32P t 1 h (/ 6 ) + 0.68P t ° h (J 6 ). 



dynamics of the qubit. A more detailed analysis of its 
damping time and period will be presented elsewhere ; 
here we focus on the amplitude of these oscillations. 

In the data shown here, the contrast is 32%, but we ob- 
served a strong dependence on the exact bias point, rang- 
ing between 5% and 40%. This visibility was lower than 
in our previous experiment [f| . We checked that the am- 
plitude lb of the current pulse was optimal by measuring 
the switching probability as a function of lb both with- 
out microwave (Ps W (^b)) and after a tt pulse (Psw(Ib))- 
Typical data are shown in figure 03. The solid black 
curve Pt h (Ib) is a fit to Ps W (h) assuming that switching 
occurs by thermal activation across the potential barrier 
0- The dotted line P? w (Ib) was calculated with the same 
parameters assuming that the qubit is always in 1. The 
curve Pg W (Ib) is then well approximated by a weighted 
sum of the curves P^ h {I) and P® h (Ib) with coefficients 
p g = 0.68 and p e — 0.32 (grey solid line). From this we 
conclude that the detection efficiency D of our detector, 
as designed, should be of at least 85%, but that the ex- 
cited state of our qubit is only populated up to 32% when 
switching to the running state occurs at the end of the 
DCP. Thus, it seems that the qubit partially relaxes to- 
wards its ground state during the measurement process. 
As we will show below, this problem does not occur with 
the alternative measurement scheme described hereafter. 

The SQUID, connected to the capacitor C s h by alu- 
minum wires of inductance L, behaves as an oscilla- 
tor with a characteristic frequency called the plasma 
frequency 2nf p = [{L + Lj)C s h]^ 1 ^ 2 ■ This frequency 
depends on the bias current and on the criti- 
cal current Ic via the Josephson inductance Lj = 
4>o / (271-/(7 \/l — I 2 Jlc) ■ Thus, the plasma frequency 
takes different values /p ' or depending on the 

qubit 's state. This effect has already been observed in 
|lOj| by measuring the transmission of a weak probe at 
a frequency close to f p while keeping the SQUID in the 
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FIG. 3: (insert) Typical resonant activation peak (width 
40 MHz), measured after a 50 ns microwave pulse. Due 
to the SQUID non-linearity, it is much sharper at low than 
at high frequencies, (figure) Center frequency of the reso- 
nant activation peak as a function of the external magnetic 
flux (squares). It follows the switching current modulation 
(dashed line). The solid line is a fit yielding the values of the 
shunt capacitor and stray inductance given in the text. 



zero-voltage state. With a probe of larger power, it has 
recently been demonstrated that a Josephson junction 
could switch between two distinct oscillation states while 
staying in the zero- voltage regime, due to its non-linearity 
. Such a bifurcation amplifier could be used to detect 
the state of a charge qubit 0. Here we apply a mi- 
crowave pulse at a frequency close to /°. We adjust the 
power so that the SQUID switches to the finite voltage 
state by resonant activation [14( if the qubit is in state |0), 
whereas it stays in the zero- voltage state if it is in state 
|1). If the two resonant activation peaks corresponding 
to the frequencies /° and /* are distinct enough, the 
switching probability of the SQUID measures the prob- 
ability that the qubit is in the state 0. Similarly, in £|, 
the quantum state of a current-biased Josephson junc- 
tion (phase-qubit) has been detected by applying a mi- 
crowave pulse that induces switching determined by the 
qubit state. In our experiment, the detector is distinct 
from the qubit. Although switching to the finite- voltage 
state is a dissipative process potentially harmful for quan- 
tum computation, the specific interest of this scheme is 
that it allows to heavily filter the bias-current line since 
no fast DC or RF pulse is needed there. 

A typical resonant activation peak is shown in the in- 
sert of figure |3J Its width depends on the frequency, 
ranging between 20 and 50 MHz. This corresponds to 
a quality factor between 50 and 150. The peak has an 
asymmetric shape, with a very sharp slope on its low- 
frequency side and a smooth high-frequency tail, due to 
the SQUID non-linearity. We could qualitatively recover 
these features by simple numerical simulations using the 
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FIG. 4: (a) Principle of the detection scheme. After the Rabi 
pulse, a microwave pulse at the plasma frequency resonantly 
enhances the escape rate. The bias current is maintained for 
500ns above the retrapping value, (b) Resonant activation 
peak for different Rabi angle. Each curve was offset by 5% 
for lisibility. The Larmor frequency was f q — 8.5 GHz. Pulse 
2 duration was 10 ns. (c) Resonant activation peak without 
(full circles) and after (open circles) a ir pulse. The continuous 
line is the difference between the two switching probabilities, 
(d) Rabi oscillation measured by DC current pulse (grey line, 
amplitude A — 40%) and by resonant activation method with 
a 5 ns RAP (black line, A = 62%), at the same Larmor 
frequency. 



RCSJ model [Tj|. The resonant activation peak can be 
unambiguously distinguished from environmental reso- 
nances by its dependence on the magnetic flux thread- 
ing the SQUID loop $ S(? . Figure shows the measured 
peak frequency for different fluxes around $ S9 = 1.5$o, 
together with the measured switching current (dashed 
line). The solid line is a numerical fit to the data using 
the above formulae. From this fit we deduce the follow- 
ing values C sh = 12 ± 2pF and L = 170 ± 20pH, close 
to the design. We are thus confident that the observed 
resonance is due to the plasma frequency. 

We then measure the effect of the qubit on the reso- 
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nant activation peak. The principle of the experiment 
is sketched in figure 0Ji. A first microwave pulse at the 
Larmor frequency induces a Rabi rotation by an angle 
0\. A second microwave pulse of duration T2 = 10ns is 
applied immediately after, at a frequency f% close to the 
plasma frequency, with a power high enough to observe 
resonant activation. In this experiment, we apply a con- 
stant bias current lb through the SQUID (lb — 2.85/iA, 
h/Ic = 0.85) and maintain it at this value 500 ns after 
the microwave pulse to keep the SQUID in the running 
state for a while after switching occurs. This allows suf- 
ficient voltage to build up across the SQUID and makes 
detection easier, similarly to the plateau used at the end 
of the DCP in the previously shown method. At the end 
of the experimental sequence, the bias current is reduced 
to zero in order to retrap the SQUID in the zero-voltage 
state. We measured the switching probability as a func- 
tion of fi for different Rabi angles 9\. The results are 
shown in figure 0Jd. After the microwave pulse, the qubit 
is in a superposition of the states and 1 with weights 
Po = cos 2 (9i/2) and p\ = sin 2 (9i/2). Correspondingly, 
the resonant activation signal is a sum of two peaks cen- 
tered at fp and with weights po and p\ , which reveal 
the Rabi oscillations. 

We show the two peaks corresponding to 9\ — (curve 
Pg W , full circles) and 9\ — ir (curve P£ w , open circles) in 
figure Hk. They are separated by f°-fp = 50 MHz 
and have a similar width of 90 MHz. This is an indi- 
cation that the tt pulse efficiently populates the excited 
state (any significant probability for the qubit to be in 
would result into broadening of the curve P£ w ), and 
is in strong contrast with the results obtained with the 
DCP method (figure^). The difference between the two 
curves S(f) = P° w — P£ w (solid line in figure ®) gives 
a lower bound of the excited state population after a ir 
pulse. Because of the above mentioned asymmetric shape 
of the resonant activation peaks, it yields larger absolute 
values on the low- than on the high-frequency side of the 
peak. Thus the plasma oscillator non-linearity increases 
the sensitivity of our measurement, which is reminiscent 



of the ideas exposed in 12]. On the data shown here, 
S(f) attains a maximum S max — 60% for a frequency f£ 
indicated by an arrow in figure 0p. The value of S max 
strongly depends on the microwave pulse duration and 
power. The optimal settings are the result of a compro- 
mise between two constraints : a long microwave pulse 
provides a better resonant activation peak separation, 
but on the other hand the pulse should be much shorter 
than the qubit damping time T±, to prevent loss of ex- 
cited state population. Under optimized conditions, we 
were able to reach S max = 68%. 

Finally, we fixed the frequency fi at the value f£ and 
measured Rabi oscillations (black curve in figure^Ji). We 
compared this curve to the one obtained with the DCP 
method in exactly the same conditions (grey curve) . The 
contrast is significantly improved, while the dephasing 
time is evidently the same. This enhancement is partly 
explained by the rapid 5 ns RAP (for the data shown in 
figure 0Jl) compared to the 30 ns DCP. But we can not 
exclude that the DCP intrinsically increases the relax- 
ation rate during its risetime. Such a process would be 
in agreement with the fact that for these bias conditions, 
Tl ~ 100 ns, three times longer than the DCP duration. 

In conclusion, we demonstrated a new method for de- 
tecting the state of a persistent-current qubit. We ex- 
ploited the dependence of the SQUID's plasma frequency 
on the qubit 's state which we measured by applying a mi- 
crowave pulse that induces resonant activation. It pro- 
vided large detection efficiency even for short microwave 
pulses (5 ns). The contrast of the Rabi oscillations was 
mainly limited by the resonant activation peaks separa- 
tion. By improving the quality factor of the plasma os- 
cillator and using composite pulses, we might thus reach 
the single-shot regime while keeping the bias current line 
heavily filtered and having a large dephasing time. 
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